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ABSTRACT 



Context. The appearance of the B[e] phenomenon in evolved massive stars such as B[e] supergiants is still a mystery. While these 
stars are generally found to have disks that are cool and dense enough for efficient molecule and dust condensation, the origin of the 
disk material is still unclear. 

Aims. We aim at studying the kinematics and origin of the disk in the eccentric binary system GG Car, whose primary component is 
|*V^ proposed to be a B[e] supergiant. 

Methods. Based on medium- and high-resolution near-infrared spectra we analyzed the CO-band emission detected from GGCar. 
The complete CO-band structure delivers information on the density and temperature of the emitting region, and the detectable 13 CO 
bands allow us to constrain the evolutionary phase. In addition, the kinematics of the CO gas can be extracted from the shape of the 
first 12 CO band head. 

Results. We find that the CO gas is located in a ring surrounding the eccentric binary system, and its kinematics agrees with Keplerian 
rotation with a velocity, projected to the line of sight, of 80 ± 1 kms~'. The CO ring has a column density of (5 ± 3) x 10 21 cm" 2 and 
a temperature of 3200 ± 500 K. In addition, the material is chemically enriched in 13 C, which agrees with the primary component 
being slightly evolved off the main sequence. We discuss two possible scenarios for the origin of the circumbinary disk: (i) non- 
conservative Roche lobe overflow, and (ii) the possibility that the progenitor of the primary component could have been a classical 
Be star. Neither can be firmly excluded, but for Roche lobe overflow to occur, a combination of stellar and orbital parameter extrema 
would be required. 



Key words. Stars: emission-line, Be - Stars: early-type - supergiants - circumstellar matter - Stars: individual: GG Car 



I/"") [ 1 . Introduction GGCar displays radial velocity and light variability in both 

* ■ the optical and UV light curves, suggestive of an eclipsing binary 

• The Galactic emission line star G GCar was discovere d mor e syste m (Her nandez et al. [19811 Gosset et al. [19841 [T9851 Brandi 

^ ; than a century ago (Pickering[l896l Pickering & Fleming! 1896J). et al. [13571 . Based on radial velocity measurements of the 

t-H • In addition to strong Balmer line emission, its optical spec- blues hifted absorptions of the Hei PCygni profiles, Marchiano 

•• ■ trum displays numerous emission lines from permitted and for- e t al. (12012b recently refined the period of this eccentric binary 

I bidden transi tions of predominantly low-ionized metals (e.g. system to 31.033 + 0.008 d. In addition, these authors obtained 

V . Swings |1974j Marchiano et al. |2012|>. Strong line profile vari- orbital parameters for the two components, as well as improved 

; ation of emission features in the o ptical (e.g. Machado et al. ste llar parameters for the primary, summarized in TableQ] The 

03 ; |2004|l and near-infrared (Morris et al. 1996) was reported, which mass G f t h e primary and its temperature estimate confirm the su- 

might be interpreted as variability in the density structure of the pe rgiant status o f GGC ar. With this knowledge, the analysis of 

wind. Furthermore, the spectral energy distribution shows a pro- Marchiano et al. (l2TjT2l also provided the distance of 5 ± 1 kpc to- 

nounced infrared excess, characteristic of emission from hot cir- war d qq Car based on two independent methods. This distance 

cumstellar dust (Allen [1973] Cohen & Barlow [1983} Bouchet estimate places GGCar in the vicinity of other evolved mas- 

& Swings [1982). These reported characteristics, together with s i ve stars such as the luminous blu e varia bles (LBVs) HRCar 

the temperature (Lopes et al. U992| l and luminosity estimates (d = 5 + I kpc; van Genderen et al. [1991] Hutsemekers & van 

(McGregor et al. 119881). resulted in the classification of GG Car Drom[l9jED and AG Car (d = 6 + 1 kpc; Humphreys et al. [19891 . 

as a B[e] supergiant (B[e]SG; Lamers et al.|1998|l, although its but keeping in mind that distance estimates generally have high 

supergiant status was long unclear due to highly uncertain dis- uncertainties. 

tance estimates. Optical polarimetric (G nedin et al. [19921 Klare & Neckel 

119771 Barbier & Swings [1982 ) as well as interferometric 



* Based on observations collected with the ESO VLT Paranal (Borges Fernandes 120 10D observations confirmed that the dust 
Observatory under programs 384.D-0613(A) and 088.D-0442(B) distribution is non-shperically symmetric, most probably located 
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Fig. 1. Fit (red/gray) to the flux-calibrated, continuum-subtracted SINFONI spectrum (black) of GGCar. Positions of the CO band 
heads and hydrogen Pfund lines are indicated. 



Table 1. Stellar parameters of the primary component of GG Car 
from Marchiano et al. (120121 1. 



T cff [K] 
L/L e 

R, [R ] 
M[M ] 

in 

d [kpc] 

E(B - V) [mag] 



23 000 ±2 000 
(2.6+ 1.0) x 10 5 
32 ±8 
26 ±4 
63 ±9 
5± 1 

0.51 ±0.15 



within a disk or disk-like structure. Preliminary analysis of the 
interferometric data deliver an inclination angle of about 60°, 
which excellently agrees with the inclination of the orbital plane 
of 54° - 72° found by Marchiano et al. d20T2l >. Furth ermore, Ha 
spectropolarimetric studies of Pereyra et al. ((2009) revealed a 
rotating disk-like structure. A disk is an ideal environment, not 
only for dust, but also for molecules, such as the very robust CO 
molecule. Emission in CO bands from GG Car has been reported 
by McGregor et al. (119881 1 and Morris et al. (119961 I. 

Kraus (2009) reported on the detection of 13 CO-band emis- 
sion from GGCar based on a low-resolution (R ^ 1500) spec- 
trum obtained with AMBER at the VLTI. The strength of the 
isotopic bands indicated that the circumstellar material is chem- 
ically enriched with processed material that must have been 
transported to the stellar surface and released into the environ- 
ment by either strong non-spherically symmetric stellar winds 
or mass ejection events. The observed enhancement with 13 C, 
however, strictly disagrees with theoretical model predictions of 
both rotating and non-rotating stellar models, and the disagree- 
ment could only be solved by the postulation of a higher mass 
(and thus luminosity) of GG Car than those reported in the lit- 
erature. The new stellar parameters determined by Marchiano 
et al. (120 12t see Table[TJ) moved the position of GGCar in the 
Hertzsprung-Russell (HR) diagram to a more appropriate place. 
Still, the proper determination of the amount of 13 CO and ac- 
cordingly the 12 C/ 13 C ratio, which is an excellent tracer of the 



stellar evolutionary phase (Kraus 2009; Liermann et al. l20101 l. is 
lacking. 

Probing the evolutionary stage of B[e]SGs is an important is- 
sue for understanding their mass-loss history, and hence the ori- 
gin and shaping mechanism of their circumstellar environment 
resulting in the B[e] phenomenon. For this, the amount of chem- 
ically enriched matter is essential, as well as the geometry and 
kinematics of the circumstellar material. The goal of the present 
paper is, therefore, to extract these based on a detailed study 
of CO-band emission in medium- and high-resolution Zf-band 
spectra of the Galactic B[e]SG star GG Car. 

2. Observation and data reduction 

A medium-resolution {R = 4500) ^-band ( 1 .95 - 2.45 fim) spec- 
trum for the Galactic B[e] star GGCar was obtained on 2012 
January 9, using the Spectrograph for INtegral Field Observation 
in the Near -Infrared (SINFONI; Eisenhauer et al. l2003l Bonnet 
et al. 120041 1 on the VLT UT4 telescope. The observations were 
performed with an AB nod pattern and an 8 x 8 arcsec 2 field of 
view. For telluric correction and flux calibration, a B-type stan- 
dard star was observed at similar airmass. The data were reduced 
with the SINFONI pipeline (version 2.2.9). 

For flux calibration, we applied a Kurucz model (Kurucz 
119931 1 to the standard star spectrum and scaled it to its Two 
Micron All-Sky Survey (2MASS, Skrutskie et al. l2006b # s -band 
magnitude. The signal-to-noise ratio (S/N) of the final spectrum 
is ~ 60. The heliocentric velocity correction was carried out and 
the spectrum was dereddened with the extinction value listed 
in TableQ] applying the interstellar reddening law of Howarth 
(HUD. 

A high-resolution spectrum (R ^ 50000) in the A'-band 
(2.276 - 2.326 yum) was obtained on 2009 December 2 using 
the CRyogenic high-resolution InfraRed Echelle Spectrograph 
(CRIRES; Kaeufel et al. 120041 1 . The observations were per- 
formed with a slit-width of 0.4 arcsec, and a standard nodding 
on-slit strategy was applied to remove sky- and detector glow. 
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Fig. 2. Fit (red/gray) to the observed, normalized CRIRES spec- 
trum of GG Car (black) considering the emission to consist of 
either pure CO band (top panel) or CO band plus Pfund line 
emission (bottom panel). 

For telluric corrections, a standard star was observed close in 
airmass and directly after the science object. The data reduc- 
tion was performed with the CRIRES pipeline. The final spec- 
trum was corrected for heliocentric velocity, and it has an S/N of 
~ 120. 

3. Results 

In addition to the emission from both the 12 CO and 13 CO bands, 
the SINFONI spectrum (Fig.Q} clearly shows emission lines 
from the hydrogen Pfund series. To model their contribution, 
we first adopted a Gaussian-like profile, indicative of an origin 
within a wind or shell. However, no reasonable fit of their profile 
shape could be achieved. 

Double-peaked lines have been identified in the optical 
spectra of GGCar, e.g., forbidden emission lines of [Oi] and 
[Can] (Muratore et al. 2012) and permitted emission lines of 
Fen (Marchiano et al. I20121 >. Furthermore, emission lines from 



the hydrogen Balmer and Paschen series appear to be double- 
peaked, typically with the red peak more intense than the blue 
one (Marchiano et al. 120121 1. Whether the Pfund lines show the 
same behavior is hard to tell because of the poor quality of the 
SINFONI spectrum (see noise level shortward of 2.29 /im). 

Pfund emission is usually formed within a narrow region. 
Considering that the line profiles could be double-peaked, the 
emission might arise from a (narrow) ring within the rotating 
disk. The spectral resolution of SINFONI, y res - 60kms~ 1 , 
smears out the possibly double-peaked profile shape so that we 
cannot measure the proper rotation velocity from the peak sep- 
aration. The rotation speed, projected to the line of sight, was 
therefore obtained from fitting the line widths with a rotation 
profile convolved to the proper spectral resolution. The optimum 
fit was achieved for u ro t,pf sin i = (94 + 2) km s _1 . 

The Pfund series reaches its edge (i.e., the Pfund discontinu- 
ity) blueward of the CO band head. Furthermore, the higher the 
transition, the shorter the wavelength separation between indi- 
vidual lines so that the lines start to blend. The broad profiles 
of the Pfund lines in our spectra together with the low spec- 
tral resolution result in the formation of a quasi-continuum of 
the higher transitions of the Pfund series, i.e., in the wavelength 
region A Z> 2.325 /mi and extending blueward of the CO band 
head (see, e.g., Kraus et al. 2000). However, the SINFONI spec- 
trum does not show any indication for such a quasi-continuum 
on the short wavelength edge of the CO spectrum, so the num- 
ber of Pfund lines in the spectrum must be limited. This agrees 
with the fact that population of the higher levels in especially 
high-density media is usually prevented by pressure ionization 
effects. The maximum Pfund transition cannot be determined 
purely based on the SINFONI data, but the highest Pfund line 
should occur redward of the CO band edge. 

The CRIRES spectrum (Fig.[2j shows the fully resolved 
structure of the first 12 CO band head emission. While fitting this 
short portion of the emission spectrum is not sensitive to the CO 
temperature and column density, it delivers detailed knowledge 
of the CO kinematics. The blue shoulder and red peak are typical 
characteristics for emission in a rotating ring or disk (e.g., Carr 
et al. [19931 Carr [T995l Najita et al. [19961 Kraus et al.[2000), and 
their separation (as shown by the arrow in the top panel of Fig.[2]i 
directly determines the disk rotational velocity, projected to the 
line of sight, in the CO-forming region. 

The observed band head characteristics could also be inter- 
preted with an equatorial outflow in the form of a gaseous ring 
expanding with constant velocity (e.g., Kraus et al. 120101 1. as 
the resulting double-peaked profile is indistinguishable from the 
rotating one. Nevertheless, we consider here the scenario of a 
rotating disk as more reliable based on the studies of Pereyra et 
al. ( 120091 . 

From our spectrum we obtain u ro t,co sin i = (80 + 1) km s . 
This value is lower than that obtained from the width of the 
Pfund lines and hence perfectly agrees with the scenario of a 
Keplerian rotating disk in which the CO-emitting region is lo- 
cated farther away than the Pfund-emitting region. From the pro- 
jected velocity of the CO gas we obtain the actual rotational ve- 
locity of y ro t,co = (91.5 + 7.5) kms -1 , assuming that the disk is 
co-planar with the orbital orientation. 

For the model computations, we made use of the CO disk 
code developed by Kraus et al. ([2000), modified to account for 
the emission of 13 CO (Kraus 2009; Oksala et al. in preparation). 
Former calculations of CO emission from a disk revealed that it 
is sufficient to consider the innermost ring of gas with a constant 
temperature, column density, and rotation velocity (e.g., Kraus 
2009 1 Liermann et al. l20T0l Cidale et al. l20T2l . 
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Table 2. Parameters of the CO disk. 



Tco [K] 
Nco [cm" 2 ] 
HI 

12 C/ 13 C 
^rot.co [km s~ 
Vmtpt [kms- 1 
Aco [AU 2 ] 



3 200 ± 300 
(5 + 3) x 10 21 
63 + 9 
15 + 5 
91.5+7.5 
108 ±11 
1.33 ±0.73 



Notes. The inclination angle is taken from Marchiano et al. ( 120 12t . 



Fitting the CRIRES spectrum with pure CO emission did not 
deliver a reasonable fit, because several features of the spectrum 
could not be satisfactorily reproduced (see top panel of Fig.|2j. 
Therefore, emission from the rotationally broadened Pfund lines 
was included. The best fit (bottom panel of Fig.[2]i was obtained 
for a cut-off in the Pfund series at the transition Pf 63. The con- 
tribution of the Pfund emission is shown in this figure as well. 

With this cut-off in the Pfund series and fixing the kinematics 
in both the CO and Pfund line emitting regions, we were finally 
able to obtain reasonable fits to the total SINFONI spectrum 
(Fig. ID- The complete set of derived CO disk parameters is listed 
in Table[2] The last row gives the area of the emitting CO ring 
considering the disk inclination and the distance of d = 5 + 1 kpc 
towards GGCar as obtained by Marchiano et al. (2012). 

The enrichment in I3 C obtained from the strength of the 
13 CO band emission of GGCar is lower than postulated by 
Kraus (2009), who suggested 12 C/ 13 C < 10. This disagreement 
can be explained by the much lower resolution of the AMBER 
spectrum in which the strongly rotationally broadened Pfund 
lines, which are even in the SINFONI spectrum fully resolved 
only at the long-wavelength end, are completely smeared out. 
Still, their presence pushes the 13 CO band heads up, falsely sug- 
gesting a stronger enrichment than the actual value of I2 C/ 13 C 
= 15 + 5 found with our detailed modeling. 

The derived CO parameters (i.e., temperature and column 
density) are similar to those found for other B[e]SGs such 
as LHA 120-S 12 and L HA120 -S73 in the Large Magellanic 
Cloud (Liermann et al. l20T0j l, L HA 1 1 5-S 65 in the Small 
Magellanic Cloud (Oksala et al. I20T21 and CPD-5 29243 
(Cidale et al. l20T2l i and MWC 349A (Kraus et al. [20001 in the 
Galaxy. Only one Galactic object, HD 327083, seems to be ex- 
treme in the sense that its disk, studied by Wheelwright et al. 
j20T2l l, has the coolest (r C o = 1721 + 540 K) and densest 
(AT C0 = (0.4 + 1.8) x 10 24 cirT 2 ) CO gas. However, in the case 
of HD 327083 only the first band head was studied, which is a 
poor tracer for density and temperature. These parameters can 
be determined with a high precision only from a complete CO 
band covering three or more band heads. 



4. Discussion 

4. 1 . Circumbinary disk of GG Car 



The binary parameters derived by Marchiano et al. (2012)) re- 
sulted in a size of the semi-major axis of the system (corrected 
for inclination) of A - 0.65 + 0.04 AU and individual stellar 
masses of M pr i m ^26 + 4 M and M sec ^12 + 2 M Q . Interpreting 
the velocity of the CO gas with rotation in a Keplerian disk, the 
distance of the CO emitting ring would be R C o = 2.88 ± 0.89 AU 
when orbiting the primary, and Rqo = 1.35 ± 0.45 AU when or- 
biting the secondary star. In both cases, this distance is larger 
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Fig. 3. Keplerian velocity distribution considering the total min- 
imum and maximum binary mass. Horizontal lines represent the 
range in rotational velocity of the CO gas. The location range of 
the CO ring is marked (dashed lines). 



than the binary separation. The CO must therefore be located in a 
circumbinary disk orbiting a central mass of Mbi na ry -38+6 M©. 

Fig. [3] displays the Keplerian velocity as a function of dis- 
tance for the minimum and maximum binary masses. The hor- 
izontal lines show the range in rotation velocity obtained from 
modeling the CRIRES data, and the intersections deliver the 
range of possible distances of the circumbinary CO-emitting ring 
from the center of gravity. From the size of the emitting CO 
area (Tabled the outer edge, R out , of the CO emitting ring can 
be calculated as a function of its inner edge, R[ n . Considering 
the error in the CO area, we computed the two limiting scenar- 
ios and obtained sizes for the CO ring that range from R out = 
2.97 +0.04 AU for the minimum inner edge of / fr^mm = 2.90 AU, 
and R out = 5.54 + 0.02 AU for the maximum inner edge of 
^in.max = 5.50 AU, indicating that the emitting CO gas is con- 
fined in a narrow ring. Considering even the largest distance, the 
CO-emitting ring falls well inside the value of 230 stellar radii 
(~ 34 AU for a primary stellar radius of 32 Rq) suggested by 
Marchiano et al. (I20121 I as the maximum possible distance of 
the CO-emitting region based on temperature arguments. 

The B[e]SG star HD 327083 was also recently found to have 
a circumbinary rotating disk displaying strong CO band emis- 
sion (Wheelwright et al. I2012I I. In both this and the GGCar 
system, the minimum inner edge of the CO ring is located at a 
similar distance of ~ 3 AU. Other interesting objects in which a 
close (in projection) binary component was discovered by means 
of interferometry are the A[e] supergiant HD 62623 (Millour et 
al. I201 II although this star was recently re-classified as a B8-9 
Ib-II star by Borges Fernandes et al. in prep.) and the B[e]SG 
candidates HD 87643 (Mill our et al. 120091 and V921 Scorpii (= 
CD-42 1 1721) (Kraus et al. |20T2T >. Only the first two have been 
found to display CO band emission, and detailed modeling will 
determine the kinematics of the CO gas and possibly help to dis- 
tinguish between the scenarios of circumstellar or circumbinary 
origin (Muratore et al. in preparation). 
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Fig. 4. HR diagram showing the position of the primary com- 
ponent of GGCar with respect to new evolutionary tracks of 
Ekstrom et al. d2012l l for stars at solar metallicity rotating ini- 
tially with Uini/ucrit = 0.4. The enrichment in 13 C over the course 
of stellar evolution is indicated. 



4.2. Evolutionary stage of GG Car 

The origin of the B[e] phenomenon in massive evolved stars like 
B [e] SGs is still not fully resolved. A disk or ring of high-density 
material containing gas and dust is usually seen around these ob- 
jects, but the disk formation mechanism is not fully understood 
and the disk kinematics have not been sufficiently studied yet. 
Because post-main sequence massive stars possess strong line- 
driven stellar winds during their main-sequence evolution, these 
disks cannot be remnants from their pre-main sequence stage, 
but likely formed from highly non-spherically symmetric stel- 
lar mass-loss. A possible scenario discussed in the literature to 
explain the formation of outflowing high-density disks in single- 
star evolution is the bistability mechanism induced by rotation 
(Lamers & Pauldrach I1991I Pelupessy et al. 2000). However, 
the density of these disks increases enormously when the slow- 
wind solution (Cure 2004) at high rotation rates (Q > 0.6) is 
considered (Cure et al. 2005). In both, the rapid stellar rota- 
tion is an important prerequisite, but it is even more important 
to know if a massive main-sequence star could remain close to 
the limit of critical rotation for a substantial fraction of its life- 
time (Langer 1998). Rotation velocities at a substantial fraction 
of their critical values have been found so far for two B[e]SGs 
LHA 115-S 23 (Kr aus et al. 120081 and LHA 1 15-S 65 (Zickgraf 
2000; Kraus et al. [20T0b in the Small Magellanic Cloud. But 
whether all B[e]SGs are indeed rapid rotators is not known. 

Other scenarios consider that binary interaction up to full 
merging of massive binaries can cause the B[e] phenomenon 
in the supergiants (Langer & Heger 1998; Podsiadlowski et al. 
2006 ). However, the binary merger scenario seems to hold so far 
for only one object of this class, i.e., the B[e]SG R4 in the Small 
Magellanic Cloud. 

The appearance of the B[e] phenomenon in the eccentric bi- 
nary system GG Car might also be linked to binary interaction, 
at least to some extent, and in the following we briefly discuss 



the evolutionary stage of the primary component of GG Car and 
its possible evolutionary history. 

Figure [4] shows the location of the primary component of 
GG Car in the HR diagram based on its effective temperature and 
luminosity (Table[TJ. Also shown are new tracks from Ekstrom et 
al. (120 121 ) for single-star evolution at solar metallicity and initial 
rotation speeds of about 40% of the critical velocity. The rota- 
tion speed might not be fully appropriate for GG Car, for which 
no rotational velocity could be obtained so far due to the lack in 
photospheric lines. But an initial rotation speed like this has been 
found as the peak of the velocity distribution for young B stars 
(Huang et al. 2010), and thus can be regarded as representative 
for the general behavior of rotating stars. The different line styles 
used along the evolutionary tracks in Fig. [4] mark regions with 
different values of the 12 C/ 13 C surface abundance ratio. This ra- 
tio starts with an initial (interstellar) value of ~ 90 on the zero- 
age main sequence, and drops to a final value of 55 5 at late, 
post-red supergiant evolutionary phases, tracing a strong enrich- 
ment in I3 C during the evolution of massive stars. The value of 
the carbon isotopic ratio found for the circumbinary material of 
GG Car is I2 C/ 13 C= 15+5. This high value excludes a late (i.e., 
post-red supergiant) evolutionary phase, but excellently agrees 
with the primary having just evolved beyond the main sequence 
(see Fig.|4]) considering that the primary has evolved unaffected 
by its companion. 

The initial mass and current location of GG Car in the HR 
diagram is very similar to another B[e]SG, the object LHA 120- 
S 12 in the Large Magellanic Cloud. Furthermore, the CO ring 
around this star was found to have very similar CO column den- 
sity (N C o = (2.5 + 0.5) x 10 21 cm" 2 ) and temperature (2800 + 
500 K), and its carbon isotopic ratio was determined to 20 ± 2 
(Liermann et al. l20T0b . Like LHA 120-S 12, GG Car is hence an- 
other example of a B[e]SG with only slight enhancement of 13 C 
and consequently is in a pre-red supergiant phase, while other 
objects were found to show strong 13 C enrichment, which clas- 
sifies them as post-red supergiants or post-yellow hypergiants 
(see e.g., Liermann et al. 120101 Muratore et al. 2010; Oksala et 
al. in preparation). 

Although the carbon isotopic ratio is the strongest indica- 
tor, there are additional arguments in favor of an evolutionary 
phase of GGCar just beyond the main sequence. Unlike the 
well-known LBVs HRCar and AG Car, no nebulosity was de- 
tected around GG Car (Thackeray 1950). Therefore, it seems less 
evolved than its close-by neighbors. Furthermore, assuming that 
the three stars HRCar, AG Car and GGCar have similar ages 
(based on similar sky location and distance), it is clear that the 
more massive stars are more evolved. Comparing the ages of 
AG Car and HR Car with the evolutionary tracks of stars in the 
proper mass range for the primary of GG Car delivers an evolu- 
tionary stage just beyond the main sequence, in perfect agree- 
ment with its position in the HR diagram. 

In summary, all evidence ( 13 C enrichment, age, missing 
LBV nebulosity) seems to speak in favor of a single-star evo- 
lutionary scenario for the primary component of GG Car, which 
is currently in an early post-main sequence evolutionary phase, 
according to the HR diagram. 

4.3. Origin of the circumbinary disk 

One question that remains to be discussed concerns the origin of 
the cool and dense circumbinary disk seen in both CO and dust 
emission. The transportation and accumulation of stellar mate- 
rial into circumbinary orbit might have happened in two different 
ways, either via non-conservative Roche lobe overflow occurring 



5 



Kraus et al.: Molecular emission from GG Car's circumbinary disk 




4.6 4.4 4.2 1 2 3 4 S 

log T„„ [K] q 



Fig. 5. Left panel: Stellar radius, normalized to the minimum bi- 
nary separation of the GG Car system, during the main and post- 
main sequence evolution of stars with different initial masses, 
plotted in terms of the effective temperature variation. Right 
panel: Relative effective Roche lobe radius for a range of mass 
ratios. The value (with errors) for GG Car is marked. 

earlier in the evolution of the primary, or the primary became a 
Be star during its late stages of main-sequence evolution. In the 
following we briefly describe and discuss both options. 

Roche lobe overflow results when the stellar radius becomes 
larger than the size of its Roche lobe. This can happen at different 
stages of stellar evolution. To compute the effective radius, r\, of 
GG Car's Roche lobe we applied the equation 

r x _ QA9g 2 ^ 

A 0.6 ? 2 / 3 +ln(l+ ? 1 / 3 ) ( ' 

provided by Eggleton (1983 ). The parameter q is the mass ra- 
tio and A is the binary separation. This equation is strictly valid 
for circular orbits only. In eccentric binaries, the mass transfer 
will take place when the two stars are closest to each other, i.e., 
during periastron passage. We can hence approximate an equiv- 
alent Roche lobe for eccentric binaries by replacing the binary 
separation A by the periastron distance A per j = A(l - e), where 
e is the eccentricity. The behavior of r\ /A pel -i as a function of the 
mass ratio q is shown in the right panel of Fig. [5] For the chosen 
range, the radius of the Roche lobe increases with increasing q. 
The mass ratio (including its errors) of GG Car and the resulting 
size of the effective Roche lobe is indicated. 

To determine if the radius of GG Car's primary star could 
have filled its Roche lobe during its past evolution, we com- 
puted the evolution of the stellar radius during and beyond the 
main sequence for stars with different appropriate initial masses. 
Assuming that the periastron distance remained constant during 
the star's evolution, we normalized the stellar radii to GG Car's 
current periastron distance and plot them versus the evolution of 
the stellar effective temperature in the left panel of Fig. [5] The lo- 
cation of GG Car's primary is included in the plot, as well as the 
range in Roche lobe radii for the current mass ratio. Obviously, 
as long as neither the binary separation nor the mass ratio has 
changed during earlier evolution, the primary of GGCar has 
not filled its Roche lobe yet. If, however, we postulate that the 



circumbinary material originates from non-conservative Roche 
lobe overflow in the past, then the original mass ratio must have 
been higher than the current value. A higher mass ratio, on the 
other hand, increases the effective radius of the Roche lobe. 
Hence, for mass transfer to happen, the original periastron dis- 
tance must have been (much) smaller than it is today, and the 
mass transfer must have resulted in widening the system. 

While this scenario cannot be excluded (see, e.g., Sepinsky et 
al. 120071 [2009), non-conservative mass transfer usually leads to 
a shrinking orbit. In addition, standard theory of binary evolution 
assumes that the originally eccentric orbit circularizes immedi- 
ately after the onset of mass transfer (e.g., Hurley et al. 2002). 
This assumption seems to be justified because only few eccentric 
evolved binary systems are observed (e.g., de Mink et al. I20071 >. 
In this context, the eccentric orbit seen for GG Car might speak 
in favor of an unevolved binary system in which no mass transfer 
has occurred yet, which would agree with the primary radius be- 
ing smaller than the effective Roche lobe size during periastron 
passage (see Fig.|5]l. 

However, Roche lobe overflow cannot be entirely excluded, 
because Fig. [5] was computed for the mean value of the perias- 
tron distance. Taking into account its error, which is based on 
the error in inclination angle and eccentricity, and considering 
the extreme case of minimum mass ratio and maximum stellar 
radius, the effective radius of the Roche lobe drops below the 
stellar radius. Therefore, Roche lobe overflow could be possi- 
ble, but would be restricted to an extreme corner in the overall 
parameter space. 

In the second possible scenario, the primary of GG Car could 
have evolved from a Be star progenitor, particularly if its mass is 
closer to the lower boundary of the primary's mass range. This 
value would agree better with those of massive Be stars found 
close to the end of and slightly beyond the main sequence. Be 
stars show Balmer emission that is supposedly arising from a 
geometrically thin, ionized gaseous disk in Keplerian rotation 
(e.g., Porter & Rivinius 2003 ). The disk is thought to be formed 
by mass loss through an equatorial, viscous decretion disk (Lee 
et al. [T99T1 Porter [T9991 Krticka et al. l20TTb triggered by th e 
close to critical rotation of the star (e.g., Townsend et al. 2004). 
Typical sizes range from a few up to a few hundred stellar radii. 

Assuming that the plane of the disk around the primary is 
coplanar with the binary orbit, this disk can cause mass loss 
into circumbinary orbits in two ways. The binary separation of 
GGCar during periastron is A peri = 0.468 + 0.048 AU, which 
corresponds to roughly 5-15 times the stellar radius during the 
main sequence evolution (see Fig. [5}. Therefore, during perias- 
tron passage even a small disk will extend beyond the Roche 
radius and disk mass could stream to the companion through 
the Lagrangian point and/or leave the Roche surface to accu- 
mulate in circumbinary orbits. On the other hand, the secondary 
will have to pass through the disk as well, and disk material ap- 
proaching the star from behind will be accelerated and hence 
kicked out of the system. 

It is hard to tell which of the two suggested scenarios is real. 
But given the need for an extreme parameter combination for 
Roche lobe overflow to work, in combination with the eccentric 
orbit of GG Car that points to an unevolved binary system, the 
Be star scenario might be slighlty more favorable. Nevertheless, 
only proper theoretical modeling of both can give a more con- 
clusive answer to the question of the evolutionary history of 
GG Car. In addition, a proper study of the kinematics of the cir- 
cumstellar/circumbinary material needs to be performed based 
on many individual tracers of different regions. The easiest indi- 
cators for kinematics are certainly the forbidden emission lines. 
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Among the observed lines are [Oi] and [Can], which are the 
most reliable tracers for the high-density disk regions (Aret et 
al. l2012l l. The profiles of the [Oi] lines, which seem to consist 
of two distinct double-peaked profiles (Muratore et al. l2012t . al- 
ready show that the distribution and kinematics of the material 
in the GG Car system is highly complex. High-quality observa- 
tions, especially over the full orbital period, are of vital impor- 
tance to extract a more complete picture of this fascinating ob- 
ject. 

5. Conclusions 

We presented our analysis of the CO-band emission of the ec- 
centric binary system GGCar. CO bands are detected in both 
medium- and high-resolution /T-band spectra obtained with the 
SINFONI and CRIRES spectrographs. The SINFONI spectrum 
extends over the full CO-band structure, which gives insight 
into the density (Nco = (5 + 3) x 10 21 cm~ 2 ) and tempera- 
ture (Tqo = 3200 + 500 K) of the emitting region, which we 
find to be confined in a narrow circumbinary ring. In addition, 
the SINFONI spectrum covers the band heads from the isotopic 
molecule 13 CO, allowing us to determine the evolutionary phase 
and hence the age of the star at the time of the material ejection. 
Our modeling derives a value 12 C/ 13 C = 15 + 5, confirming a 
slightly evolved (supergiant) nature of GG Car's primary com- 
ponent. The CRIRES spectrum contains the spectrally resolved 
structure of the first 12 CO band head. From its shape we find that 
the motion of the CO gas is consistent with Keplerian rotation 
with a velocity, projected to the line of sight, of 80 + 1 km s _1 . 

Based on these results, we discussed the origin of the cir- 
cumbinary disk in this eccentric binary system. In particular, 
two scenarios were taken into consideration. The first one con- 
sidered the possibility that the primary component underwent 
non-conservative Roche lobe overflow during its past evolution, 
while the second considered that the mass transport and accu- 
mulation on circumbinary orbits could have happened during Be 
phases of the primary component. We cannot firmly prove or 
disprove the validity of either scenario. However, the conditions 
for non-conservative Roche lobe overflow are only met in an 
extreme corner of the parameter space. Therefore, slight pref- 
erence is given to the suggestion that during late phases of its 
main-sequence evolution, GG Car's primary could have been a 
classical Be star. 
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